Electronic Excited States in Bilayer Graphene Double Quantum Dots 
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We report tunneling spectroscopy experiments on a bilayer graphene double quantum dot device 
that can be tuned by all-graphene lateral gates. The diameter of the two quantum dots are around 
50 nm and the constrictions acting as tunneling barriers are 30 nm in width. The double quantum 
dot features addition energies on the order of 20 meV. Charge stability diagrams allow us to study 
the tunable interdot coupling energy as well as the spectrum of the electronic excited states on a 
number of individual triple points over a large energy range. The obtained constant level spacing 
of 1.75 meV over a wide energy range is in good agreement with the expected single-particle energy 
spacing in bilayer graphene quantum dots. Finally, we investigate the evolution of the electronic 
excited states in a parallel magnetic field. 
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Graphene quantum dots (QDs) are interesting candi- 
dates for spin qubits with long coherence times [1]. The 
suppressed hyperfine interaction and weak spin-orbit cou- 
pling [2I, Q make graphene and flat carbon structures in 
general ^promising for future quantum information tech- 
nology [j] . Significant progress has been made recently in 
the fabrication and understanding of graphene quantum 
devices. A "paper-cutting" technique enables the fabrica- 
tion of graphene nanoribbons (5l-[l3|, quantum dots [3- 
[l9| , and double quantum dot devices pQM23j, where a 
disorder-induced energy gap allows confinement of indi- 
vidual carriers in graphene. These devices allowed the 
experimental investigation of excited state s |l6l . l2l| , spin 
states [m and the electron-hole crossover [l^ However, 
all of these studies were based on single-layer graphene 
and showed a number of device limitations related to 
the presence of disorder, vibrational excitations and to 
the fact that the missing band gap makes it difficult to 
realize soft confinement potentials and "well-behaving" 
tunneling barriers. In particular, it has been shown that 
intrinsic ripples and corrugations in single-layer graphene 
can lead to unintended vibrational degrees of freedom (2^ 
and to a coherent electron-vibron coupling in graphene 
QDs [ill. Bilayer graphene is a promising candidate 
to overcome some of these limitations. In particular it 
allows to open a band gap by an out-of-plane electric 
field [26|42S|, which may enable a soft confinement po- 
tential and may reduce the influence of localized edge 
states. More importantly, it has been shown that rip- 
ples and substrate-induced disorder are reduced in bilayer 
graphene [29|, which increases the mechanical stability 
and suppresses unwanted vibrational modes. 

Here, we present a bilayer graphene double quantum 
dot (DQD) device, with a number of lateral gates. These 
local gates allow to tune transport from hole to electron 
dominated regimes and they enable to access different de- 
vice configurations. We focus on the DQD configuration 
and show characteristic honeycomb-like charge stability 
diagrams, with addition energies on the order of 20 meV. 
Most interestingly, we observe electronic excited states, 
which are equally spaced over a wide energy range, in 
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FIG. 1: (color online) (a) Scanning force microscope image 
of the investigated double quantum dot device. The scale 
bar is 200 nm. (b) 2D-line of the Raman spectrum fitted 
with four Lorentzians. The separation of the inner peaks 
measures 18.8 cm~^. (c) Back gate characteristics recorded 
at Vh = 20 mV and T = 1.3 K. The arrow marks the BG 
regime where the following measurements have been taken, 
(d) Schematic illustration of the effective band structure of 
the device highlighting the three tunneling barriers (hatched 
areas) induced by the local constrictions (see panel a). 



good agreement with the expected single-particle con- 
finement energy in bilayer graphene QDs. 

The DQD device shown in Figure la is fabricated based 
on bilayer graphene obtained from mechanical exfolia- 
tion of bulk graphite. The bilayer graphene flakes are 
deposited on a highly p-doped silicon substrate with a 
295 nm silicon oxide layer. Raman spectroscopy mea- 
surements are used to unambiguously identify bilayer 
graphene. In Figure lb we show the 2D Raman peak 
recorded on the flake which has been used to fabricate 
the device shown in Figure la. In contrast to single-layer 
graphene we find that four Lorentzians are required for 
fitting the 2D line shape. In particular, the left shoul- 
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der (see arrow in Fig. lb) and the spacing between 
the inner two peaks, which is 18.8 cm~^ (see dashed 
hnes in Fig. lb) provide a clear fingerprint for bilayer 
graphene [30|, l3l| . Electron beam (e-beam) lithography 
is used to pattern the device etch mask on the isolated 
graphene fiake. Reactive ion etching based on Ar/02 
plasma is employed to transfer the pattern to the bi- 
layer graphene. Finally, an additional e-beam and lift-off 
step is used to contact the nanostructure with Cr/Au 
(5 nm/50 nm) electrodes. 

Figure la shows a scanning force micrograph of the 
investigated bilayer DQD. The diameters of the etched 
quantum dots (QDs) measures roughly 50 nm while the 
width of the 100 nm long constrictions leading to the 
dots measure 30 nm. For tuning the individual tunnel- 
ing barriers and QDs, the device includes lateral bilayer 
graphene gates positioned in a distance less than 80 nm 
next to the active structure. In this study we used the 
left and right gate (LG and RG) to change the number 
of carriers in the left and right QD (LQD and RQD), 
respectively, while the central gate (CG) and the second 
left gate (LG*) are used to tune the inter-dot coupling 
and the left barrier, respectively. Additionally, the back 
gate (BG) is used to adjust the overall Fermi level. 

The measurements have been performed in a pumped 
^He system with a base temperature of 1.3 K and in a di- 
lution refrigerator with an electron temperature around 
100 mK. We have measured the two- terminal conduc- 
tance through the bilayer graphene DQD by applying a 
symmetric dc bias voltage Vb while measuring the current 
through the device with a resolution better than 50 fA. 

The source-drain current of the device measured at 
finite bias (Vb = 20 mV) over a large BG voltage range 
is shown in Figure Ic. In close analogy to single-layer 
graphene nanodevices we observe a region of suppressed 
current separating the hole (left inset) from the electron 
transport regime (right inset). This so-called transport 
gap [8[, which extends for the investigated device from 
roughly 32 V to 48 V is expected to be mainly caused by 
the local tunneling barriers formed by the three 30 nm 
narrow constrictions (see illustration in Figure Id) [3, 
flT^ . The reproducible sharp conductance resonances in 
and around the transport gap region are due to localized 
states in the constrictions whereas the overall hole-doping 
is most likely due to atmospheric O2 binding (32| . 

In order to access the DQD regime we fix the BG volt- 
age to a value inside the transport gap, such that the 
source and drain potentials are within the valence band 
(see arrow in Figure Ic and Fermi level depicted in Fig- 
ure Id). In Figure 2a we show a measurement of the cur- 
rent as function of the central gate voltage and the left 
and right gate voltages (Vig* = Vrg) while keeping the 
BG voltage fixed (arrow in Fig. Ic). Two characteristic 
slopes are observed which separate areas of suppressed 
current from the area of elevated current in the lower 
left corner. This measurement suggests that the energy 
diagram shown in Figure Id (including a dominant cen- 
tral barrier) is a useful description for the three different 
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FIG. 2: (color online) (a) Source-drain current as a function 
of the central gate voltage (Vcg) and the voltages applied to 
the outer side gates (Vz^* ^Vrg) at Vb — 20 mV and T — 1.3 K. 
(b,c) Schematic band structure depicting the pure hole trans- 
port regime (2) and a single quantum dot regime with lifted 
central barrier (3), respectively. The dashed lines denote the 
position of the Fermi level, (d) Charge stability diagram at 
finite bias {Vh = 10 mV) recorded in the double quantum 
dot regime (regime (1), see also Figure Id) as highlighted by 
the black bar in panel (a), (e) Local resonances in the right 
constriction as function of the left and right gate voltages 
measured in regime (3) with additionally lifted left tunnel- 
ing barrier (V^^* = and Vig < 0), such that only the right 
tunneling barrier crosses the Fermi level {Vb — 10 mV) |33| . 

transport regimes (l)-(3) found in the data. In regime 
(2) all three side gates are tuned to very negative volt- 
ages, such that all three effective tunneling barriers are 
pushed in the conduction band and hole transport takes 
place throughout the whole structure (c.f. Fig. 2b). By 
increasing the voltages on the left and right gate we en- 
ter regime (3), which leads to a significantly decreased 
current. In this configuration the Fermi level crosses the 
two tunneling barriers induced by the two outer constric- 
tions (c.f. Fig. 2c). Thus the device is in a single quan- 
tum dot regime. Finally, if we increase the voltage on 
the central gate as well, the current is almost completely 
suppressed [region (1)]. The Fermi level crosses in this 
configuration all three effective energy gaps and the de- 
vice is in the DQD regime (c.f. Fig. Id). The LQD and 
RQD are consequently defined by the three energy gaps 
as shown in Figure Id. The two slopes separating the dif- 
ferent regimes are determined by the relative lever arms 
of the lateral gates (CG, LG* and RG) on the three tun- 
neling barriers. While the two outer tunneling barriers 
can be controlled well by the left, right and the central 
gate {a = AVig* ^rg/ ^Vcg ^ 0.55), the dominating cen- 
tral barrier can only be weakly tuned by Vig* and Vrg 
{j3 = /SVcg/ /SVig* ^rg ^ 0.13), which is in good agreement 
with the geometry of the device. 

By fixing the central gate voltage we can record charge 
stability diagrams in different regimes by varying the left 
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and right gate voltages independently. In Figure 2d we 
show a finite bias (14 = 10 mV, T = 1.3 K) measure- 
ment taken in regime (1) (see black bar in Fig. 2a). A 
honeycomb pattern characteristic for the charge stability 
diagram of a DQD can be observed [34|. Electrical trans- 
port through the DQD is only possible in the case where 
the energy levels in both dots are aligned within the 
source-drain bias window. In close analogy to single-layer 
graphene, the overall transmission through the DQD is 
modulated due to local resonances in the narrow con- 
strictions. This can be seen best in Figure 2e, where 
we show a charge stability diagram recorded in regime 
(3) (c.f. Fig. 2c) but with lifted left barrier {Vig* = 
and Vig < 0). Thus the Fermi level only probes localized 
states in the right tunneling barrier. In this measurement 
we observe only one dominating slope {AVig/AVrg ^ 0.2, 
dashed line), which moreover corresponds to the slope of 
the transmission modulation in Figure 2d [33| . In regions 
of high transmission (e.g. lower left corner in Fig. 2d) the 
connecting lines between the triple points become visible. 
Along these lines, only one of the dot levels is within the 
bias window, leading to current by cotunneling processes. 
The relative lever arms between the left and right gates 
acting on the two dots (LQD and RQD) are determined 
from the slopes of these cotunneling lines (see dashed 
lines in Fig. 2d) delimiting the hexagons, = 0.85 

and (^['g ,^g = 0.12. Please note that in this system the 
RQD can almost equally be tuned with the KG as with 
the LG, which again is consistent with the device geom- 
etry. 

A high-resolution close-up of Figure 2d, but with re- 
versed bias, Vb = -10 mV is shown in Figure 3a. The 
bias-dependent extensions of the triangular-shaped re- 
gions allows the determination of the conversion factors 
between gate voltages and energies (see illustration in 
Fig. 3b) [33 • The lever arm between the left (right) 
gate and the left (right) dot is = V^/SVig = 0.056 
{a^g = 0.019) 0,133. This allows to extract the single- 
dot addition energies = af^ • AVig = 23 meV and 
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try in the QD island sizes. The two quantum dots are 
considered to be mainly formed by the width-modulated 
bilayer graphene structure and the local disorder poten- 
tial. In particular, the disorder potential may have sig- 
nificant impact infiuencing the extend and location of the 
individual tunneling barriers |8|]. In Figure 3c we show 
a similar measurement from a different cool down. The 
BG voltage has been fixed in the center of the trans- 
port gap. We observe slightly modified dot sizes with 
addition energies of E^= 18 meV and E^= 21 meV. 
Interestingly, the extracted addition energies are in rea- 
sonable agreement with values from single-layer graphene 
quantum dots with a similar size [19]. In the following 
we focus on this more symmetric regime, which moreover 
has been investigated in a dilution refrigerator with an 
electron temperature of around 100 mK. The increased 
energy resolution makes the edges of the triangles much 



FIG. 3: (color online) (a) High-resolution close-up of Fig- 
ure 2d but with reversed bias (T^^-IO mV; T=1.3 K) high- 
lighting the double dot characteristic honeycomb pattern, (b) 
Schematic charge stability diagram denoting all quantities 
necessary to deduce the gate lever arms, addition energies and 
the mutual capacitive coupling energy (for more information 
see text), (c) Similar measurement as in panel (a) but at a 
different cool down and an electron temperature of around 
100 mK. The well defined triple points allow to extract the 
inter-dot coupling energy (see arrows) and additional excited 
states can be observed. 



more defined, such that the mutual capacitive coupling 
energy (£^^) between the two QDs can be extracted 
quantitatively. The coupling energy between both dots 
can be determined from the splitting of the triple points, 
E"^ = afg ■ AVl^ = a^g ■ AV^g=2A meV. (see Fig. 3b 
and arrows in Fig. 3c). 

By changing the central gate voltage we are able to 
change the couphng energy between the left and right 
QD. Figures 4a and 4b show close-ups of the very same 
triple point (14 = 6 mV) for two different inter-dot cou- 
pling energies E'^^ 4.8 meV (^^=7.42 V, Fig. 4a) and 



. AV^g = 13 meV, which refiect an asymme- £:^= 3.0 meV (14^=7.54 V, Fig. 4b), where the dif- 



ferent triple point separations can be seen. Figure 4c 
presents a more detailed analysis of how the coupling en- 
ergy changes with small variation in V^g for three individ- 
ual nearby triple points (the examples shown in Figures 
4a and 4b are highlighted by arrows). While we observe 
an increase of the coupling energy for the neighboring 
triple points [see (N,M) and (N+1,M), where N and M 
are the number of carriers on the left and right dot, re- 
spectively] for increasing Vcg we observe a decrease of 
the coupling energy by more than a factor of 2 for the 
triple point (N+5,M-1). This is in agreement with ear- 
lier studies on graphene QDs [20|, [111 showing a strongly 
non-monotonic dependence of the conductance on gate 
voltages due to the sharp resonances in the constrictions. 
Interestingly, the extracted couphng energies are roughly 
by a factor of 2 larger, as compared to earlier single-layer 
DQD studies [lO,!!!!, which might be related to the con- 
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FIG. 4: (color online) (a,b) Very same triple points in a weaker (Vcg = 7.54 V, Vb = 6 mV) and stronger (Vcg — 7.42 V) 
inter-dot coupling regime respectively, (c) Mutual capacitive coupling energy depending on the central gate voltage Vcg, 
recorded for three individual triple points. N, M denote the number of electrons occupying the dots [see also inset in panel 
(d)]. (d) High-resolution close-up of Fig. 3c. The triple points show a dense set of electronic excited states. Inset: Schematic 
energy diagram showing the ground states of the double system and the first excited state of the left dot. (e,f) Different triple 
points featuring similar excited state spectra {Vb = 10 mV and 6 mV respectively), (g) Single- level spacing as a function of 
the relative change of the carrier occupation number on the left dot. The circles show the experimental data and the dashed 
horizontal line marks the experimental mean value (A=1.82 meV). The three curves represent the model for a single- layer QD 
for different absolute occupation numbers (N = A^o + AN) while the solid horizontal line marks constant level spacing for a 
bilayer dot. The QD diameter is assumed to be 50 nm. 



siderably smaller size of the investigated device. 

The increased energy resolution also reveals additional 
fine structure inside the individual triple points. This, 
for example, can be seen in Figure 3c. The data ex- 
hibit distinct lines of increased conductance parallel to 
the base line of the triple points. Figure 4d shows a 
close-up of Figure 3c highlighting a set of electronic ex- 
cited state resonances parallel to the triple points base 
line (see arrows in Fig. 4d). Along such a parallel line the 
inter-dot detuning energy, £, is constant and the current 
is increased due to transport through the excited state 
(see dashed lines in the inset in Fig. 4d). The electronic 
nature of these excited states can be shown by the mag- 
netic field dependence as discussed below. The detuning 
energies of these 5 visible excited state resonances are 
found to be £ = 1.7, 3.3, 4.9, 6.5, 8.1 meV, respectively 
(see black arrows in Fig. 4d). Figures 4e and 4f present 
high-resolution scans of different triple points recorded 
at different bias voltages. The triple points in Figure 4e 
show excited states parallel to the base line as well as a 
resonance parallel to its right edge. The first two excited 
states are again found at detuning energies of e = 1.7 and 
3.8 meV (see black arrows). The conduction suppression 
near the right edge can be attributed to variations of the 
coupling between energy levels of the left dot and the 
source lead Q. In Figure 4f we show a triple point pair 
recorded at Vb = -6 mV with very similar excited state 



energies {e = 1.8 and 3.5 meV; see black arrows). 

In total we analyzed more than 50 excited states in a 
wide energy range. The electronic excited state energy 
spacing has been found to be constant over the entire 
range and a value of A = 1.75 ±0.27 meV has been 
extracted. Interestingly this value is in good agreement 
with the single-particle confinement energy in disk-like 
bilayer graphene QDs, which can be estimated by us- 
ing the density of states for bilayer graphene D{E) = 
7i/7r(fi'UF)^, where vp is the Fermi velocity and 71 = 
0.39 meV is the inter-layer hopping energy. Consequently 
the single-particle level spacing is given by 



where d is the diameter of the bilayer QD. Assuming a 
QD diameter of d = 50 nm, which is in reasonable agree- 
ment with the lithographically defined QD we obtain a 
constant (carrier number independent) level spacing of A 
= 1.71 meV, which is in good agreement with the exper- 
imental data. 

The charge carrier number independent level spacing is 
in contrast to single-layer graphene QDs where the single- 
particle level spacing, AfTV) = ?iVF/d\fN^ depends on 
the number of carriers, N [16|- In Figure 4g we show the 
single-level spacing for a 50 nm diameter single-layer and 
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bilayer graphene QD as function of the excess number 
of carriers on the left dot, AN = N - Nq. The sohd 
horizontal line marks the constant level spacing for the 
bilayer QD, whereas the three l/>/7V-dependent curves 
belong to a single-layer QD with different absolute num- 
bers of charge carriers on the dot {Nq = 5, 25 and 50). 
The experimental data, which are depicted by circles in 
Fig. 4g have been extracted from charge stability dia- 
grams as shown e.g. in Fig. 3c. For larger excess carrier 
numbers (AA^ > 40) we also made use of relative gate 
lever arms to estimate AN. Since we cannot fully ex- 
clude charge rearragements, these data have to be taken 
with care. Please note that the total carrier numbers (A/", 
M) on the two quantum dots are unknown. In Fig. 4g we 
plot only excited states energies extracted from measure- 
ments at Vb < and we focus on changes of the number 
of carriers on the left QD only (A7V, see also inset in 
Fig. 4d). For this limited set of excited states we find a 
constant value of A = 1.82 ±0.14 meV, which is depicted 
by the dashed horizontal line in Fig. 4g. Please note 
that a constant level spacing might also result from an 
effective constant density of states induced by disorder. 
However, following Ref. [17], we estimate the effective 
dot disorder potential to be on the order of the addition 
energy of 5-10 charge carriers. Since we are probing a 
significantly larger energy range we can rule out disor- 
der to be cause of the observed constant level spacing, 
and attribute it to true band structure properties. Con- 
sequently, the model suitable for describing single-layer 
QDs [A (A/") oc 1/a/]V] does not fit the experimental data 
for any reasonable Nq whereas the model for a bilayer 
QD [A(A^) = const.] is in good agreement. This leads 
to the conclusion that the investigated QDs are indeed 
extending over both graphene sheets forming the bilayer 
system. 

Figure 5 shows a study of the evolution of excited 
states as function of the magnetic field oriented paral- 
lel to the bilayer graphene plane proving the electronic 
nature of these states. In Figures 5a and 5b we show a 
pair of triple points (Vb = -6 mV) for two different val- 
ues of magnetic field. Excited states parallel to the base 
line as well as one parallel to the left edge of the triangle 
are visible (dotted lines). For zero magnetic field again a 
level spacing of 1.8 meV (see e.g. arrow in Fig. 5a) can 
be observed. Figure 5b shows the very same triple points 
at a parallel magnetic field {B\\) of 2 T. The position of 
the triangles in the gate voltage plane almost does not 
change at all (see crosses and solid lines in Fig. 5a, b). 
This holds even for -fields up to 9 T (as shown in 
Fig. 5c) and shows that the effects on the orbital parts of 
the wave functions are unaffected. The parallel magnetic 
field has two effects on the triple points, i.e. triangles. 
First, the conductance at large detuning values slightly 
increases with increasing -field. This is particularly 
true for the left triangle, where even parts of the triangle 
(left edge) are strongly suppressed at low 5|| -fields. Since 
the associated slope differs from the actual triple point 
edge we attribute this effects to a modulated transmis- 
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FIG. 5: (color online) (a,b) The very same individual triple 
points recorded ai B\\ = T and B\\ = 2 T respectively 
(Vb — —6 mV). (c) Current as a function of the detuning 
energy [measured along the dashed line indicated in panel 
(b)] for different 5|| -fields up to 9 T. Traces are offset by 
2 pA for clarity, (d) Current as a function of the detuning 
energy and the 5|| -field recorded at a different triple point. 
The corresponding triple point is shown in the inset, where 
the scale bars are 50 mV in gate voltages. 



sion from the drain to the right dot {Vb < 0). Second, 
the 1*^ excited state parallel to the base line splits into 
two separated peaks (see dotted lines in Fig. 5b). This 
peak splitting is linear in magnetic field up to roughly 
5 T as shown in Figure 5c, where we show current line 
cuts along the detuning energy axis (see dashed line in 
Fig. 5b). The different traces are recorded for different 
parallel magnetic fields (see labels at the right side) and 
they are offset by 2 pA in current for clarity. We ob- 
serve that the position of the base line (^=0) and the 
2nd excited state {e = 3.6 meV) are constant up to 9 T 
as function of the -field (see vertical dashed lines in 
Fig. 5c). In contrast, the 1^^ excited state splits into two 
peaks, where the characteristic linear peak splitting is 
0.2 meV/T (see diverging dashed lines). A similar peak 
splitting can also be observed on different triple points. 
In Figure 5d we show for example measurements taken on 
a triple point, where the 3^^ excited state {e = 5.3 meV) 
is splitting. Here the current at Vb = -8 mV is plotted 
as a function of detuning energy and parallel magnetic 
field. The inset shows the corresponding triple points as 
well as the detuning axis (dashed line). While the posi- 
tion of the baseline and the 1^^ excited state are constant 
in 5|| -field, the 2^^ excited state shows a down-shift in 
energy and the 3^^ excited state clearly splits linearly in 
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two peaks (see dashed lines). The base hne as weh as the 
1** and 2^^ excited state completely vanish at increased 
parallel magnetic fields. The observed B-field dependent 
peak shifts have again a slope of either or ± 0.1 meV/T. 

All observed peak splittings are likely to be originated 
from Zeeman spin splitting. This is motivated by the lin- 
ear 5|| -field dependence of the peak splittings (up to high 
fields) and their characteristic energy scale, which is in 
reasonable agreement with g/j^s = 0.116 meV/T, assum- 
ing a ^-factor of 2. The ^-factor 2, moreover, would be 
in agreement with recent spin-resolved quantum interfer- 
ence measurements in graphene [35| and observations of 
spin states in a single-layer graphene quantum dot [l9|. 
However, a detailed analysis of the different splitting and 
non-splitting (e-dependent) inter-dot transitions is diffi- 
cult since the carrier number, including the total spin 
and the degree of lifting the valley degeneracy are not 
known. Possible transitions leading to peak splittings 
may be traced back to increased total spin numbers, non- 
spin conserving inter-dot tunneling or different g-factors 
in both quantum dots. Moreover, we can not exclude that 
the local ^-factor is increased by an unintentional layer- 
dependent doping or by the local gate potentials, which 
both could potentially lead to an out of plane electric 



field, which is expected to increase the effective ^-factor 
in bilayer graphene [36]. 

In summary, we have fabricated and characterized a 
tunable bilayer graphene double quantum dot based on 
a width-modulated graphene nanostructure with lateral 
graphene gates. Its functionality was demonstrated by 
the observation of a tunable inter-dot coupling energy 
and the electronic excited state spectra over a wide en- 
ergy range. We have shown that - in contrast to single- 
layer graphene - the single-particle level spacing is inde- 
pendent of the number of charge carriers on the bilayer 
quantum dots. By applying a 5|| -field we observed inter- 
dot transition energy splittings on the order of Zeeman 
splittings. These results give insights into tunable bilayer 
graphene double quantum dot devices and open the way 
to study individual spin states and spin coherence times 
in more detail in future experiments. 
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